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Abstract 

t — \ The nonperturbative nature of nucleon-nucleon interactions as a function of a mo- 

mentum cutoff is studied using Weinberg eigenvalues as a diagnostic. This investi- 
ng . gation extends an earlier study of the perturbative convergence of the Born series 
to partial waves beyond the 3 Si- 3 Di channel and to positive energies. As the cut- 
off is lowered using renormalization-group or model-space techniques, the evolution 
of nonperturbative features at large cutoffs from strong short-range repulsion and 
the iterated tensor interaction are monitored via the complex Weinberg eigenvalues. 
When all eigenvalues lie within the unit circle, the expansion of the scattering ampli- 
tude in terms of the interaction is perturbative, with the magnitude of the largest 
eigenvalue setting the rate of convergence. Major decreases in the magnitudes of 
*^ ' repulsive eigenvalues are observed as the Argonne v\s, CD-Bonn or Nijmegen po- 
tentials are evolved to low momentum, even though two-body observables are un- 
changed. For chiral EFT potentials, running the cutoff lower tames the impact of 
the tensor force and of new nonperturbative features entering at N 3 LO. The effi- 
cacy of separable approximations to nuclear interactions derived from the Weinberg 
^ analysis is studied as a function of cutoff, and the connection to inverse scattering 
is demonstrated. 
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1 Introduction 



Perturbative expansions of the Lippmann-Schwinger equation for the scatter- 
ing amplitude using conventional nucleon-nucleon interactions fail to converge 
in low partial waves for several reasons. First is the strongly repulsive short- 
range interaction present in all conventional potential models, which implies 
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large contributions from poorly constrained high-momentum components in 
intermediate states and requires a nonperturbative treatment. Second is the 
tensor force, e.g., from pion exchange, which is singular at short distances and 
requires iteration in the triplet channels. Finally, the presence of low-energy 
bound states or nearly bound states in the S- waves imply poles in the T matrix 
that render the perturbative Born series divergent or nearly so. These features 
greatly complicate the nuclear few- and many-body problems. 

However, the nucleon-nucleon potential is not an observable, fixed uniquely 
from experiment, and these nonperturbative features can vary widely among 
the infinite number of potentials capable of accurately describing low-energy 
physics. The effective field theory (EFT) philosophy leads us to exploit this 
freedom by choosing inter-nucleon interactions appropriate for the physics 
problem of interest. In this regard, a compelling alternative to conventional 
potentials for the many-body problem is to use energy-independent nucleon- 
nucleon interactions with a variable momentum cutoff, which are evolved to 
low momentum using renormalizat ion-group or model-space methods, together 
with consistent cutoff-dependent many-nucleon interactions [1,2,3,4,5,6]. These 
low-momentum potentials, generically called "Mow it") preserve the on-shell T 
matrix for momenta below the cutoff. 1 

Changing the momentum cutoff shifts contributions between nuclear interac- 
tions and loop integrals over intermediate states such that two-nucleon observ- 
ables are unchanged. This affects the perturbativeness of an interaction and 
the strength of the associated correlations in the wave functions, which means 
that many-body calculations can be simplified, particularly since Pauli block- 
ing eliminates at moderate densities (well below nuclear matter saturation) 
the impact of the bound and nearly bound states in the S-waves. In Ref. [5], 
it was shown that low-momentum interactions make nuclear matter calcula- 
tions perturbative in the particle-particle channel and that the corresponding 
low-momentum three-nucleon interactions drive saturation. This refutes the 
long-standing conventional wisdom that the nuclear matter problem must be 
nonperturbative in the inter-nucleon interactions in order to obtain saturation. 

In order to monitor and quantify our observations about perturbative expan- 
sions, a method introduced by Weinberg [9] that focuses on eigenvalues of the 
operator Gq(z)V proves very useful. This method was applied in Ref. [5] at the 
deuteron binding energy to the S-waves in free space and in the medium. In 
the present work, we provide a complete analysis in free space. The Weinberg 
eigenvalues obtained for nucleon-nucleon interactions are studied as a func- 
tion of the cutoff for various partial- wave channels and for a range of energies. 



In this work, we use V\ ow k potentials with sharp cutoffs on the relative momenta. 
The use of smooth cutoffs for V\ ow k, which have technical advantages for some 
applications, is discussed in Refs. [7,8]. 
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When these complex eigenvalues lie within the unit circle, the T matrix ex- 
pansion in terms of the potential V is perturbative, and the largest eigenvalue 
determines the rate of convergence. By focusing separately on the "attractive" 
and "repulsive" eigenvalues (see Sect. 2), we can isolate the different sources 
of nonperturbative behavior in nuclear forces. 

In Sect. 2, we apply this analysis to Vi ow fc derived from the Argonne fig po- 
tential [10], which has been used in the most accurate ab-initio calculations 
of nuclei and nuclear matter to date. The removal of the strongly repulsive 
core with decreasing cutoff is tracked quantitatively via the largest repulsive 
Weinberg eigenvalues. Comparisons to other potentials (CD-Bonn [11] and 
Nijmegen [12]) shows the model dependence inherent in the short- distance 
descriptions, which is removed as the cutoff is lowered to 2 fm~ . Chiral EFT 
potentials [13,14], which provide a more systematic and model-independent 
approach to inter-nucleon interactions, typically have cutoffs below 3.5 fm -1 , 
but running the cutoff lower tames the impact of the tensor force and of new 
nonperturbative features entering at N 3 LO. The eigenvectors corresponding 
to the Weinberg eigenvalues lead naturally to separable approximations to 
the nucleon-nucleon interaction [15]. In Sect. 3, we show that the expansion 
in separable potentials becomes more effective for lower cutoffs, and compare 
^lowfc with separable inverse-scattering solutions. Sect. 4 summarizes our con- 
clusions. 



2 Weinberg Eigenvalue Analysis for Low-Momentum Interactions 

The eigenvalue analysis introduced by Weinberg [9] provides quantitative con- 
ditions for the perturbative convergence of the T matrix. This allows us to 
identify momentum cutoffs at which repulsive core scattering and iterated 
tensor contributions become perturbative. Here we briefly review the Wein- 
berg formalism and refer the reader to the literature [9,15,16,17] for more 
details. 

The Lippmann-Schwinger equation for the T matrix with energy-independent 
potential V is given in operator form by 



T(z) = V + VG (z)T(z), 



(1) 



where the free propagator G (z) with kinetic energy H is 



G {z) 



1 



(2) 



z-H, 



o 



The T matrix can be expanded in the perturbative Born series as 



T(z) = V + VG (z)V + VG (z)VG (z)V + ... . 



(3) 
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We are interested in the convergence of this series at different energies, which 
is closely related to the presence of bound states or resonances. Suppose that 
V supports a bound state at z — E b < 0. Since T(z) has a pole at z = E b 
and the individual terms in the Born series are all finite, this implies that the 
series diverges as z — > E b . 

We can rewrite the eigenvalue problem for the bound state \b), 

(H + V)\b) = E b \b), (4) 

as 

^—V\b) = G (E b )V\b) = \b) , (5) 

which motivates the generalization to a new eigenvalue problem for complex 
energies z: 

G {z)V\^ v {z))=r lv {z)\^ v {z)). (6) 

Here the index v labels the Weinberg eigenvalues (which are always dis- 
crete [9,17]) and eigenvectors. The r\ v {z) are analytic functions cut along 
^ z < oo [17]. For positive energies, the eigenvalues are defined as the 
limits on the upper or lower edges of the cut: 

]imG (E±ie)V\*?\E)) = vl ±] (E) \¥±\E)} . (7) 

e— »U 

In the following, it is understood that for positive energy r) v (E) = r)^(E). 

As demonstrated in Ref. [9], the perturbative Born series for T(z) diverges 
if and only if there is an eigenvalue with \r] u (z)\ ^ I. The necessity of this 
condition is clear, since 

T(z) \* v (z)) = (l + Vu (z) + Vu (z) 2 + ...)V \* v (z)) (8) 

diverges for |^(-z)| ^ 1. Any scattering state \^/) can be written as a linear 
combination of the Weinberg states \^ u (z)), thus the Born series will converge 
only if all f] u (z) lie within the unit circle. Furthermore, the rate of convergence 
is controlled by the largest |77„(z)|, with smaller magnitudes implying faster 
convergence. 

A rearrangement of Eq. (6) gives a simple interpretation of the eigenvalues 
r\ v {z) in terms of the Schrodinger equation, 

(H + -ir7)|^( Z )) = z|^)>. (9) 

T lv\ z ) 

The eigenvalue r] u (z) can thus be viewed as an energy-dependent coupling that 
must divide V to produce a solution to the Schrodinger equation at energy z. 
If V supports a bound state at z — E B , then there is some v with r) u (E B ) = 1, 
which implies a divergence of the Born series for nearby energies. However, 
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what matters for convergence at a given energy z is not simply the presence of 
nearby physical bound states, but rather the entire set of eigenstates that can 
be shifted to z when the interaction is divided by i] u (z). For negative energies, 
a purely attractive V gives positive rj u (z) values, while a purely repulsive V 
gives negative eigenvalues, as the sign of the interaction must be flipped to 
support a bound state. For this reason, we follow convention and refer to 
negative eigenvalues as repulsive and positive ones as attractive. In the case 
of conventional nuclear interactions, the repulsive core generates at least one 
large and negative eigenvalue that causes the Born series to diverge in the low 
partial waves. 

For energies greater than zero, the eigenvalues are still discrete but become 
complex. Weinberg also showed that there are at most a finite number of 
eigenvalues with magnitudes greater than unity [9]. This is used in Weinberg's 
"quasiparticle" method to systematically isolate the nonperturbative parts 
of the potential in separable form. The eigenvalues vary continuously with 
energy, so they can be plotted as trajectories in the complex plane, as in Figs. 1 
and 2. For each energy, the overall largest eigenvalue determines whether or 
not the Born series converges at that energy. But by considering attractive and 
repulsive eigenvalues separately we can isolate the contributions of different 
sources of nonperturbative behavior. 

Here we calculate the Weinberg eigenvalues in different partial waves as a func- 
tion of energy for Vj wfc over a wide range of cutoffs. The construction of V\ ow k 
using renormalization-group (RG) or model-space methods is described in de- 
tail elsewhere [2,3]. In short, the RG approach starts with the partial- wave 
Lippman-Schwinger equation for the T matrix with a cutoff A on interme- 
diate relative momenta. By imposing the condition dT{k' , k; k 2 )/dA = on 
the half-on-shell T matrix, we ensure that low-momentum observables are in- 
dependent of the scale A and generate energy-independent potentials. This 
leads to coupled differential equations for the interaction V\ ow k{k',k), which 
are solved on a momentum grid suitable for Gaussian integration, starting 
from a large initial A appropriate for the starting "bare" interaction. This 
starting point could be the Argonne Vi$, CD-Bonn, or Nijmegen potential, the 
chiral EFT potentials at N 3 LO, or any other potential that reproduces two- 
body observables. For sharp cutoffs A ^ 2fm _1 , the resulting V\ ow k preserves 
the elastic phase shifts up to laboratory energies = 2A 2 ^ 330 MeV (in 
units where h 2 /m = 1). V\ ow k(k',k) is obtained on a momentum grid, where 
the eigenvalue problem Eq. (6) is given by 

2 j k2df£ Vlo, k (k',k) = ^ {k ^ u{p ^ . (10 ) 

^ 

We solve Eq. (10) by converting it to a left-eigenvalue problem, where we 
integrate over the singularity using l/(x + ie) = V/x — m5{x). Alternatively, 
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Fig. 1. Trajectories of the largest repulsive Weinberg eigenvalue in the 1 So channel 
as a function of energy for Vj ow /c derived from the Argonne v±s potential. Our 
results for selected cutoffs are indicated by the different symbols. The positions of 
the symbols on each trajectory mark the eigenvalues for center-of-mass energies 
E cm = k 2 = 0,25,66,100 and 150 MeV, starting from the filled symbol at OMeV. 
The trajectory with stars are eigenvalues for the N 3 LO potential of Ref. [13]. 
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Fig. 2. Trajectories of the largest repulsive Weinberg eigenvalue in the 3 Si- 3 Di 
coupled channel as a function of energy for V| wfc derived from the Argonne v±s 
potential. Our results for selected cutoffs are indicated by the different symbols. The 
positions of the symbols on each trajectory mark the eigenvalues for center-of-mass 
energies E cm = 0, 25, 66, 100 and 150 MeV, starting from the filled symbol at MeV. 
The trajectory with stars are eigenvalues for the N 3 LO potential of Ref. [13]. 
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Fig. 3. Trajectories of the largest repulsive Weinberg eigenvalue in the 1 So channel 
and the 3 Si- 3 Di coupled channel as a function of energy for Vi ow fc derived from the 
CD-Bonn (left) and Nijmegen I potentials (right). Our results for selected cutoffs 
are indicated by the different symbols (see Fig. 1 for the legend). The positions 
of the symbols on each trajectory mark the eigenvalues for center-of-mass energies 
E cm = 0, 25, 66, 100 and 150 MeV, starting from the filled symbol at OMeV. 



we can solve the complex right-eigenvalue problem for VGq(z) [V\^ u (z))] = 
r\ v {z) ty\fy v {z))\, which has the same spectrum as Go(z)V and integrates over 
the singularity directly. 



In Figs. 1 and 2, we show the trajectories of the largest repulsive eigenvalue 
for the S-waves as a function of positive energy over a range of cutoffs from 
A = 10fm _1 to 2fm~ 1 . The eigenvalues are for the V\ ow k evolved from the 
Argonne v is potential, and for A > 10 fm~ there are essentially no changes in 
the eigenvalues. In all cases, a trajectory starts on the real axis with energy 
E — (for repulsive eigenvalues ^(0) < 0) and evolves counter-clockwise [9]. 
If E were decreased to negative values, the trajectory would continue along 
the real axis with decreasing magnitude of r] as E becomes more negative. 
It is evident that the magnitude of the largest repulsive eigenvalue at all 
energies is greatly cutoff dependent. The decrease with lowering the cutoff 
reflects the elimination of the repulsive core of the potential. In the singlet 
channel, Fig. 1, the trajectory lies completely inside the shaded unit circle for 
cutoffs near 4fm _1 and below, which implies that the Born series becomes 
perturbative with respect to the repulsive core (but still converges very slowly 
at low energies due to the nearly bound state at threshold, see the discussion 
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on the attractive Weinberg eigenvalue and Fig. 8 below). By A = 2 far 1 the 
largest repulsive eigenvalue is very small at all energies, which shows that the 
nonperturbative behavior from the repulsive core has been eliminated and has 
little impact on the convergence. 

In the triplet channel, Fig. 2, where the tensor contribution is active, the repul- 
sive core is still the dominant source for the largest repulsive eigenvalue. When 
the cutoff is lowered to eliminate the core, the repulsive tensor interaction is 
left as the now-dominant feature in the triplet channel, so that the largest re- 
pulsive eigenvalues are just within the unit circle for A = 3 fa -1 and still have 
a magnitude of rj m 0.5 for A = 2 far 1 . However, the decrease from A = 3 far 1 
to 2 far 1 is significant for ensuring a convergence of particle-particle ladders 
in the nuclear medium, since second-order tensor contributions are known to 
strongly excite intermediate-state momenta peaked at k ~ 2.5 — 3.5 far in 
nuclear matter [18]. In Fig. 3 we show the largest repulsive Weinberg eigen- 
values of Viowfc obtained from the CD-Bonn and Nijmegen I potentials. The 
eigenvalues for these energies are nearly identical for cutoffs below A = 3 far 1 , 
whereas the coalescence of all momentum-space matrix elements of Vj ow fc oc- 
curs for lower cutoffs A ~ 2 far 1 [1,2]. We further observe that the singlet 
channel eigenvalues are very similar even for larger cutoffs, but the triplet 
eigenvalues differ substantially for large cutoff. 

Chiral EFT interactions typically have cutoffs below 3.5 far 1 and therefore 
are expected to be soft potentials. However, as shown in Fig. 1, the N 3 LO 
potential of Entem and Machleidt [13] has substantial repulsive eigenvalues 
in the singlet channel, which are close to the Vi ow fe results for A = 4 far 1 . 
Therefore, we study the behavior of chiral EFT interactions in more detail in 
Figs. 4 and 5, using Weinberg eigenvalues as a diagnostic. Our results show 
that the repulsive eigenvalues are large for the N 3 LO potential of Epelbaum 
et al. [14], and that there is a major decrease as the cutoff is run down from 
A = 3 far 1 to 2 far 1 . This decrease is important, as the magnitude of the 
largest repulsive eigenvalue increases in the 1 So channel for E > (see Fig. 1) 
and because the tensor force in the 3 Si- 3 D! coupled channel is softened when 
chiral interactions are evolved to lower cutoffs (see also Fig. 2). 

In addition to the tensor contribution at all orders, there are new nonpertur- 
bative sources due to the central parts of pion exchanges that enter at N 3 LO. 
This can be seen from the S-wave component of the deuteron wave function at 
N 3 LO, Fig. 13 in Ref. [14], where the very short-range part for A = 550 MeV 
is similar to what is found with a repulsive core (even for A = 450 MeV, the 
S-wave deuteron wave function lies between the bare Argonne v is potential 
and Viowfc for A = 4 far 1 ). This observation and the comparison in Fig. 4 
to the eigenvalues of the N 2 LO potential [19] suggest that the large repul- 
sive Weinberg eigenvalues at N 3 LO are due to nonperturbative features in 
the central part of the sub-sub-leading 27r-exchange interaction. This might 
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Fig. 4. The largest repulsive Weinberg eigenvalues for E = in the 1 So channel (left) 
and the 3 Si- 3 Di coupled channel (right) as a function of cutoff for V\ ow k derived 
from chiral interactions. Results are shown for the N 3 LO potential of Entem and 
Machleidt [13], for the N 3 LO potential of Epclb aum et al. [14] with different cutoffs 
A/ A (as indicated in MeV), and for the N 2 LO potential [19]. For comparison, we 
have plotted the largest repulsive Weinberg eigenvalues for Vj ow & derived from the 
Argonne vi$ potential. 
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Fig. 5. The largest repulsive Weinberg eigenvalues for E = in the 1 So channel (left) 
and the 3 Si- 3 Di coupled channel (right) as a function of cutoff for V\ ow k derived 
from chiral interactions. Results are shown for the N 3 LO potential of Entem and 
Machleidt [13] and for the N 3 LO potential of Epelb aum et al. [14] with different 
cutoffs A/ A (as indicated in MeV). 
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be understood from the short-distance behavior of the isoscalar central 2tt- 
exchange two- loop contributions, discussed in [20]. The repulsive eigenvalues 
of the N 3 LO potential of Entem and Machleidt [13], which uses A = 500 MeV, 
are smaller and close to the A = 450 MeV N 3 LO potentials of Epelbaum et 
al. [14]. The differences in implementation between [13] and [14], such as the 
fitting of low-energy constants and the regularization of two-pion exchange at 
N 3 LO are reflected in the eigenvalues. As seen in Fig. 5, however, the patterns 
of eigenvalues are not completely systematic and require further study, which 
is deferred to a future publication. In any case, Figs. 4 and 5 clearly show that 
it is advantageous to evolve chiral EFT interactions to lower cutoffs using the 
RG, since the singular central (at N 3 LO) and the singular tensor parts are 
softened for lower cutoffs. 

In Figs. 6 and 7, the magnitudes of the largest repulsive eigenvalues for the 
Argonne i>i 8 potential are shown for a set of representative partial waves. 
The results are plotted as a function of energy for selected cutoffs and as a 
function of cutoff for selected energies, respectively (where not plotted, the 
eigenvalue was too small in magnitude to unambiguously identify). The same 
monotonic decrease in magnitude with decreasing cutoff is seen in every chan- 
nel. The crossover from the nonperturbative (\r]\ > 1) to the perturbative 
regime < 1) happens at different cutoffs in different channels, but always 
for A > 3fm~ x . For completeness, we note that the second-largest repulsive 
eigenvalue is below unit magnitude in all channels for all cutoffs, but shows 
the same rapid decrease to very small values by A = 2 fm _1 . 

Next, we study the largest attractive Weinberg eigenvalues. Their trajectories 
for the S- waves are plotted as a function of energy in Figs. 8 and 9. The values 
r)(E — 0) ~ 1 reflect the presence of a bound state close to zero energy; if 
the potential were scaled in each case to V/r)(0), there would be a bound 
state exactly at zero energy. Conversely, all nuclear interactions have a unit 
eigenvalue r]{E^) = 1 in the triplet channel, corresponding to the deuteron 
binding energy = —2.225 MeV. We find that the attractive eigenvalues 
increase insubstantially in magnitude as A decreases. This small increase is 
similar for all energies. As expected, the eigenvalue close to unit magnitude 
persists as we lower the cutoff. This is because it is dictated by physics, namely 
the near-zero-energy bound states. However, we have shown in Ref. [5] that 
the large attractive eigenvalues decrease rapidly as a function of density due 
to Pauli blocking. In contrast, the large repulsive eigenvalues for large cutoffs 
are little affected by Pauli blocking. 

The role of the tensor force can be seen by comparing the 3 P and 3 Px chan- 
nels in Figs. 7 and 10. The comparison shows similar repulsive eigenvalues (in 
magnitude) at large cutoff, where the dominant effect is the short-range repul- 
sion. At the smallest cutoffs, 3 Po has tiny repulsive eigenvalues while the 3 Pi 
repulsive eigenvalues are non-negligible. This is consistent with tensor contri- 
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Fig. 6. The magnitude of the largest repulsive Weinberg eigenvalue as a function of 
center-of-mass energy E cm in selected channels for a range of cutoffs. The symbols 
label the same cutoffs and the N 3 LO potential [13] as in Fig. 1. 
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cutoff A in selected channels for three center-of-mass energies. 
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Fig. 8. Trajectories of the largest attractive Weinberg eigenvalue in the 1 So channel 
as a function of energy for V\ ow k derived from the Argonne v±s potential. Our 
results for selected cutoffs are indicated by the different symbols. The positions of 
the symbols on each trajectory mark the eigenvalues for center-of-mass energies 
E cm = 0, 25, 66, 100 and 150 MeV, starting from the filled symbol at OMeV. 
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Fig. 9. Trajectories of the largest attractive Weinberg eigenvalue in the 3 Si- 3 Di 
coupled channel as a function of energy for V| wfc derived from the Argonne v±s 
potential. Our results for selected cutoffs are indicated by the different symbols. The 
positions of the symbols on each trajectory mark the eigenvalues for center-of-mass 
energies E cm = 0, 25, 66, 100 and 150 MeV, starting from the filled symbol at MeV. 
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Fig. 10. The magnitude of the largest attractive Weinberg eigenvalue as a function 
of A in selected channels for three laboratory energies. There are no significant 
attractive Weinberg eigenvalues for any A in the Pi, 3 Pi and 3 F3 channels and 
where not plotted, the eigenvalue was too small in magnitude to unambiguously 
identify. 

butions being attractive in 3 P and repulsive in 3 Pi. Figure 10 shows that the 
largest attractive eigenvalues for channels besides the S-waves are where the 
long-range one-pion exchange tensor force is attractive ( 3 Po, 3 P2~ 3 F2, 3 D2 and 
3 G4). The eigenvalues in the channels not shown are very small in magnitude. 
The behavior of the attractive eigenvalues is similar for the other potential 
models. 



3 Separable Approximations to the Nucleon-Nucleon Interaction 



The Weinberg eigenvalue analysis leads naturally to separable approximations 
for nucleon-nucieon potentials, which may be of practical interest for solving 
few-body problems such as four-body scattering [21,22]. The Weinberg eigen- 
states {j^i/)} can be used as a basis for a separable representation of the 
potential [9,15,23,24], with various possible ways to construct such a repre- 
sentation. Here we do not attempt to optimize such an expansion, but focus 
on the efficacy of a given separable expansion as a function of the cutoff, and 
choose for illustration the expansion in Weinberg eigenvectors at zero energy. 
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For E ^ 0, we use the following separable representation 

oo 

V=Y I V\* V (E)){* V (E)\V, (11) 

where we choose the normalization of the Weinberg states such that 

(p v ,{E)\V\V v {E)) = 5 v , v . (12) 

This separable representation is easily verified by substituting it into the eigen- 
value problem, Eq. (6). The T matrix is then given by 

~ V\* V {E)){* V {E)\V fl3 x 

As in Refs. [9,15], Eqs. (6) and (13) lead to a sum rule for the exact T matrix 
at energy E, 



Tr 



G (E) T(E) G (E) T(E)} = ]T - 2 . (14) 

(l - l/Vu(E)) 



The summation in the corresponding sum rule to a rank-n approximation 
runs only from v — 1 to v — n. For E ^ each term in the sum is positive def- 
inite, so the sum rule is most saturated and the difference between and T 
minimized by choosing the n eigenvalues in order of ascending (l — l/r] u (E)j . 
This illustrates a difference between the convergence of the separable expan- 
sion and that of the Born series (which depends on the eigenvalues of largest 
magnitude rather than those closest to unity). 

For practical calculations, we can apply the rank-n approximation at fixed 
energy E to calculate the T matrix at all energies E' . In this case, the analysis 
of convergence is no longer simply given by the eigenvalues, but we expect the 
general trends to carry over. Here we choose E — and study the accuracy 
of observables calculated from the eigenvectors corresponding to the n largest 
eigenvalues, using the representation Eq. (11). As the rank n increases, all 
eigenvalues that are not included decrease, and thus our expansion improves. 
In addition, we examine how the rank-n approximation performs as a function 
of the cutoff, where V\ ow k is evolved from the Argonne v±s potential for all 
results presented in this section. 

In Figs. 11 and 12, the convergence of our truncation is shown by the relative 
error in the calculated phase shifts for E^ = 100 MeV and 250 MeV (which 
typify the full range of energies). Our results demonstrate that using a rank-n 
separable approximation at the particular energy E — MeV does reproduce 
phase shifts reasonably well for E' ^ E. Except in a few cases, the convergence 
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Fig. 11. Relative error in the phase shifts for Elab = 100 MeV as a function of A 
based on separable potentials generated from the Weinberg eigenvectors at E = 
corresponding to the largest n = 3, 5 and 10 eigenvalues. 
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Fig. 12. Relative error in the phase shifts for E^ = 250 MeV as a function of A 
based on separable potentials generated from the Weinberg eigenvectors at E = 
corresponding to the largest n = 3, 5 and 10 eigenvalues. 
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improves monotonically as more terms are included in the separable approxi- 
mation and, for a particular rank-n approximation, lower cutoffs yield better 
convergence. As expected, we observe that the convergence is slower in low 
partial waves where the tensor force is active ( 3 Si and 3 D 1 in Figs. 11 and 12). 
Equation (14) implies that the fall-off of Weinberg eigenvalues should be tied 
to the rate of convergence. This is consistent with Fig. 13, which shows this 
fall-off in the % channel for £ lab = 100 MeV. 

Further insight into how the convergence varies with the cutoff comes from 
the predicted triton binding energy. Since we are interested in convergence 
rather than absolute predictions, we consider for simplicity and convenience 
the expectation value of the Hamiltonian at each A in a basis of harmonic 
oscillators up to iV max = 40 with fixed oscillator parameter b = 1.7 fm (that 
is, the energy is not minimized with respect to this parameter). The result 
for the triton binding energy E t with the full potential is compared to results 
with different rank separable approximations at the same cutoff. The absolute 
errors are plotted in Fig. 14 as a function of rank for a variety of cutoffs. (Note 
that full predictions for the triton will vary with the cutoff, because we do not 
include the A-dependent consistent three-body force, but for our purposes this 
is not relevant.) 

At very low rank, the errors are not systematic, but by rank n = 4 there is 
a systematic decrease in the error with decreasing cutoff. Our results can be 
compared to those from the separable expansions of low-momentum potentials 
considered in Ref. [25], which used polynomial expansions up to a specified 
rank. They found as well that calculations of the triton converge to the exact 
result at lower rank as the cutoff is decreased, with convergence to four digits 
or better at rank 10, 14 and 20 for A = 3, 5 and 10 fm -1 , respectively. Overall, 
the dramatic improvement in convergence of separable expansions when cutoffs 
are lowered motivates a future study of how to optimize such expansions. This 
study should also consider alternative formulations, such as that of Ernst, 
Shakin and Thaler, which when applied to conventional (and therefore high 
cutoff) potentials at rank 4, leads to deviations of less than 50keV for the 
triton binding energy [26]. 

In addition to the separable expansion in terms of Weinberg eigenvectors, one 
can consider a separable inverse scattering approach [15,27] that is based on 
the scattering phase shifts and does not require a potential model. Here we 
make a simple exploratory study, using the separable inverse scattering so- 
lution for a single channel to compare the inverse scattering low-momentum 
interaction V inv to V\ ow k in the S and P-waves. We therefore neglect the cou- 
pling of the 3 Si and 3 P2 channels to the 3 Di and 3 F2 waves respectively. The 
generalization to coupled channels and to higher rank separable interactions 
(such as an attractive and repulsive contribution) is straightforward [27]. For 
low cutoffs, the perturbative regime with Weinberg eigenvalues \rj\ < 1 sug- 
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Fig. 14. Absolute error of the triton binding energy for separable expansions of dif- 
ferent rank. The results for different cutoffs are obtained in a fixed basis of harmonic 
oscillators (based only on two-body interactions). The error is with respect to the 
exact result for the same A and model space N max . 
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Fig. 15. Diagonal matrix elements V\ ow k(k,k) compared to separable inverse scat- 
tering results V\ nv (k,k) for A = 2.0 fm -1 in S and P-waves. 

gests a direct connection to the phase shifts, except for the S-waves with 
near-zero-energy bound states. Our results show that the inverse scattering 
solution may provide a good starting point to make this connection for sepa- 
rable interactions. 

For a single channel, the solution to the inverse scattering problem with a one- 
term separable interaction V\ nv {k, k') = gv(k) v(k') is given by (where the sign 



of g determines the sign of the interaction and v(k) = JV inv (k, k)/g) [15,27] 



Vmv{k, k) 



'k 2 + k 2 \ sin 5(k) 
k 2 ) k 



-A(fc) 



(15) 



The factor (k 2 + k^)/k 2 is present for channels with a bound state (in the 3 Si 



(16) 



partial wave fcb = ym \ E&\/h 2 with deuteron binding energy E d ) and 

2 



A(Jfe) = -V I 

71 J 



P 



Pdp xf . 



Here we have generalized the conventional inverse scattering solution (where 
the phase shifts over all energies are known) to an effective theory where all 
phase shifts are known up to the cutoff. Therefore V^ nv (fc, k) is defined only 
for momenta k < A. 
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In Fig. 15, we compare the resulting inverse scattering interaction Vi nv (k, k) 
with V\ ow k(k, k) for a representative low-momentum cutoff A = 2.0fm _1 in 
the S and P-waves. We find that the separable inverse scattering solution 
reproduces qualitatively all features of V\ owk . The largest differences are ob- 
served in the S-waves, which is due to the rank-1 approximation. Note that 
the (coupled channel) 3 P2 partial wave is well described in our single-channel 
approximation, since the coupling to the 3 F2 wave is weak. In addition, the 
inverse scattering solution leads to the same characteristic dependences close 
to the cutoff, which makes it clear that these are sharp cutoff artifacts. (Note 
that the separable inverse scattering solution would go smoothly to zero if one 
were to modify the phase shifts by a smooth cutoff regulator.) Finally, it is 
interesting that the single-channel inverse scattering solution leads to a good 
description of V\ ow k in channels where the interaction changes from attractive 
to repulsive, as is the case in the 3 P partial wave. For low-momentum in- 
teractions, where the phase shifts are known over the range of validity of the 
effective theory, the inverse scattering approach therefore presents a promising 
addition to separable expansions in Weinberg eigenvectors. 



4 Conclusions 

In this paper, Weinberg eigenvalues are used as a diagnostic to study the non- 
perturbative nature of nucleon-nucleon interactions as a function of a momen- 
tum cutoff. As the cutoff is lowered using renormalizat ion-group or model- 
space techniques, the evolution of nonperturbative features at large cutoffs 
from strong short-range repulsion and iterated tensor interactions are moni- 
tored via the repulsive complex Weinberg eigenvalues. Major decreases in the 
magnitudes of repulsive eigenvalues are observed for all conventional poten- 
tials when evolved to low momentum, even though two-body observables are 
unchanged. 

When all eigenvalues for a given channel lie within the unit circle, the T ma- 
trix expansion in terms of the potential V is perturbative in that channel, 
with the magnitude of the largest eigenvalues setting the rate of convergence. 
The repulsive eigenvalues satisfy this criterium in all channels for cutoffs of 
2fm~ , where the crossover from the nonperturbative (|^| > 1) to the per- 
turbative regime (|?7| < 1) happens for A > 3fm _1 . The largest attractive 
eigenvalues, however, remain close to or beyond the unit circle in the S-waves, 
because they are due to the bound or nearly bound states. Consequently, 
this is the only surviving source of nonperturbative behavior for the Born 
series with low-momentum interactions. As shown in Ref. [5], the large at- 
tractive eigenvalues decrease rapidly due to Pauli blocking as the density is 
increased toward saturation density. This opens the possibility of perturbative 
and therefore systematic calculations of nuclear matter. 
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Chiral EFT interactions are also low-momentum interactions and are expected 
to be soft. However, at N 3 LO the isoscalar central part of the sub-sub-leading 
27r-exchange interaction is strongly repulsive at short distances [20]. This leads 
to large repulsive Weinberg eigenvalues for N 3 LO potentials in both singlet and 
triplet channels. As a result of these nonperturbative features at N 3 LO, the 
S-wave component of the deuteron wave function [14] has a short-range part 
that is similar to what is found with a repulsive core. In addition, the tensor 
force, e.g., from pion exchange, is singular at short distances. We have found 
major decreases of the repulsive Weinberg eigenvalues as the cutoff is run 
down from A = 3fm _1 to 2fm _1 , which tames the impact of the tensor force 
and of the nonperturbative features at N 3 LO. This motives the evolution of 
chiral interactions to lower cutoffs to improve their effectiveness for few- and 
many-body problems [5,8]. 

The decrease of the Weinberg eigenvalues implies that separable expansions of 
nuclear interactions become more effective for lower cutoffs. This was demon- 
strated for NN scattering and the triton binding energy: for cutoffs of 2fm~ 1 , 
the typical relative error in the NN phase shifts is 0.1%, and the convergence 
to the exact triton binding energy is to four digits at rank 6, 10, 12 and 14 for 
A = 2,3,4 and 5fm _1 , respectively. Such separable expansions may be useful 
for four-body scattering [21,22] or to develop microscopic interactions that 
can be handled by the current methods in nuclear reaction theory. Perturba- 
tive Weinberg eigenvalues also suggest a closer connection of low-momentum 
interactions to phase shifts. This was explored with rank-1 separable inverse 
scattering solutions, which were found to qualitatively reproduce the momen- 
tum dependences of V\ ow k- 
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